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The Bloch spinors, energy spectrum and spin density in energy bands are studied for the two-dimensional 
electron gas (2DEG) with Rashba spin-orbit (SO) interaction subject to one-dimensional (ID) periodic elec- 
trostatic potential of a lateral superlattice. The space symmetry of the Bloch spinors with spin parity is 
studied. It is shown that the Bloch spinors at fixed quasimomentum describe the standing spin waves with 
the wavelength equal to the superlattice period. The spin projections in these states have the components 
both parallel and transverse to the 2DEG plane. The anticrossing of the energy dispersion curves due to the 
interplay between the SO and periodic terms is observed, leading to the spin flip. The relation between the 
spin parity and the interband optical selection rules is discussed, and the effect of magnetization of the SO 
superlattice in the presence of external electric field is predicted. 

PACS: 73.21.Cd, 85.75.-d 



INTRODUCTION 

In the past years, an increasing attention has been 
drawn to the spin related phenomena in semiconductor 
structures. This research area has developed in the new 
branch of condensed matter physics and spin electron- 
ics. The problem of spin-dependent quantum states and 
'transport phenomena in these systems are currently at- 
Itracting a lot of interest also due to their potential for 
'future electronic device applications. 

In two-dimensional semiconductor heterostructures 
the spin-orbit interaction is usually dominated by the 
Rashba coupling 1 stemming from the structure in- 
version asymmetry of confining potential. The low- 
dimensional semiconductor structures with SO interac- 
tion were studied theoretically in numerous papers [2 - 
S], including the ID periodic systems with SO couphng 

MM 

The effects of spin splitting in 2DEG were investi- 
.gated experimentally by the magnetotransport studies, 
in particular, by Shubnikov - de Haas oscillations [H]. 
The role of other SO terms, such as Dresselhaus term, 
can be estimated, for example, by optical methods 0. 
As it was demonstrated by the experiments, in the SO 
structures interesting effects may be observed, such as 
the spin Hall effect ^U] and the spin-galvanic effect 

In the present paper we study quantum states and 
the electron spin distribution in a system combining 
the spin-splitting phenomena caused by the SO inter- 
action and the external gate-controlled periodic electric 



potential. We thus want to investigate the spin ori- 
entation and spin polarization that can be achieved in 
currently manufactured gated semiconductor structures 
with lateral surface superlattice. For example, the ID 
superlattice can be fabricated by the metal gate evapo- 
ration with typical period of 50 — 200 nm. We use the 
value of lateral period in the x-direction to be a = 60 
nm which gives us the energy scale n'^fi^ /2ma^ of the 
order of 2 meV for the effective mass m* = 0.067rno 
in GaAs. The values of Rashba coupling constant for 
the most important semiconductors are in the range of 
(1...5)-10~^^ eVm. It is known also that the Rashba 
coupling strength can be modified by the gate field by 
up to 50% 12 . Below in our calculations we use the 
value a = 5 • 10^^^ eVm which gives the typical shift of 
the parabolic dispersion curves kso to be of the order 
of TT / a. So, in the structure studied in the manuscript 
the electron kinetic energy i:'^h^k'^ /2ra will be compara- 
ble to the Rashba energy ak which makes the effects of 
SO interaction and periodic potential distinguishable. 
It should be mentioned also that the energy scale stud- 
ied in our paper means that the effects discussed in the 
manuscript can be clearly observed experimentally at 
helium temperatures. 



1. QUANTUM STATES 

The Hamiltonian of our problem 



^'e-mail: khomitsky@phys.unn.ru 



H = H. 
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is the sum of the Rashba Hamiltonian Hq with the SO 
couphng strength a, 



P 

2m' 



(2) 



and the one-dimensional periodic potential V{x) of a ID 
superlattice with the period a. We choose the simplest 
form of the periodic potential 



27rcc 

V{x) = Vo cos , 

a 



(3) 



where the sign and the magnitude of Vb can be con- 
trolled, for example, by an external gate, leaving us a 
wide interval of possible Vq amplitudes with an order of 
several meV. 



11 Perturbation approach 

Qualitatively the formation of the SO-split bands 
can be seen in the perturbation approach applied in the 
problem of quantum states in quantum wires ^ and 
in the tight-binding approximation for the SO super- 
lattices 13 . Namely, a.t ky = the Hamiltonian is 
written as + V{x) — adyPx- Considering the SO 
terms —aayP^ as a perturbation, one can choose the 
zero-order wavefunction as 



±i 



(4) 



where ipmkix) are the eigenstates of the Hamiltonian 
2^ -t- V{x) and m is the band number. Thus, in the 
first order of the perturbation theory the energy spec- 
trum will be determined by the expression 



Em{kx, ky — 0) 



, , hkso dSm 
s™(M±^^ (5) 



where em{kx) is the band spectrum in the ID periodic 
potential V{x) without the SO interaction. Here we 
have used the relation :^{'ipmk \ Px \ i^mk) = ^ for 
derivation of Eq[Sl An example of the energy band spec- 
trum of two lowest bands at /cy = each double-split by 
the SO interaction is given in Fig^ One can see that 
since the derivative vanishes at the border of the 



E(ir,,ir ) [arb.u.] 




FigH Energy bands at fcy = double-split by the 
SO perturbation. The degeneracy is not lifted at ~ 
0, ±7r/a, ky = since the derivative of the un- 
perturbed band dispersion vanishes at these points. 



Brillouin zone (BZ) kx = ±7r/a. 



BZ center, the degeneracy is not lifted in these points, 
being lifted at finite ky only. 



12 Bloch spinors 

In the presence of both SO coupling and periodic poten- 
tial we construct the two-component eigenstate of the 
Hamiltonian 



^lix) 
y^lix) 



(6) 



as a superposition of two-component spinors which are 
the eigenstates of the Rashba Hamiltonian The 
wavevectors of the basis states in this superposition are 
shifted by the reciprocal lattice vector b of the supelat- 
tice: 

27r 



k„ = k -I- nb = { kx-\ n, 

a 



(7) 



n = 0,±1,±2, 
the form 



The eigenstate in band m thus has 



V2 \ Ae*' 



, A = ±1 (8) 



where k is the quasimomentum in the ID Brillouin zone 
and On = arg[fcy — iknx]- After substituting the wave- 
function JHJ into the Schrodinger equation the coeffi- 
cients a™„ are determined by the standard eigenvalue 
problem 



E)S 



AnA' 



AA' 

n' I ^ rin' 



V 



X'n' 



= 0, (9) 



and at the where E? is the energy of a free Rashba quantum state 



(10) 



Spin-orbit lateral superlattices: energy bands and spin polarization in 2DEG 



3 



taken at the point in k-space defined by {T)), i.e. E^-)^ = 
Exi^x + 2TTn/a, ky). The matrix elements in the system 
© have the form 

F„5 = K,A„„,(1 + AAV('^"-'^"')), ^^^^ 

Ann' — \Sn,n'±lj 71 = 7l' ± 1. 

The structure of matrix elements determines the 
classification of energy bands and gaps in the SO su- 
perlattice. The dependence of the matrix elements Hll() 
on the quantum numbers and ky can be obtained 
directly from Ea. (|ll(l . The matrix elements t^^,7±i 
and V~^^i describe the coupling between the nearest- 
neighboring states O with n' = n ± 1 and with the 
opposite indices A ^ A' labeling different up- and down- 
Rashba bands (|1U|I . Here the superscript ± labels the in- 
dex A = ±1. We shall see below that these elements pro- 
duce the energy gaps located inside the Brillouin zone. 
The elements V^^j^-^ and l^„7n±i ^'^^ responsible for the 
coupling between the states of the same Rashba index 
A = A'. Such elements could open gaps on the borders 
of the Brillouin zone k^ — ±7r/a. However, one can 
see from that they vanish for the case of pure elec- 
trostatic periodic potential. Then, for the values of at 
ky >• 7r/a the elements with opposite indices A ^ A', i.e. 
V^~j^i and V~^^i decrease to zero while the elements 
with A = A' approach their maximum values. 

13 Energy spectrum 

The energy spectrum calculated with the help of sys- 
tem @ is shown in Fig[21 Here in Fig[2Ji. we give an 
example of energy spectrum E{kx) at fixed ky = Q and 
in FigEh the fcj,-dependence of the same spectrum for 
kx = 0. In accordance with the Kramers theorem the 
symmetry Em]i}^) = Emx(— k) together with the sym- 
metry kx^y — *■ —kx,y takes place, and we thus show the 
spectrum only at positive kx and ky. One can see that 
the spin degeneracy at fcj, = is not lifted at the center 
and at the borders of the BZ kx = ±7r/a. This result 
was also obtained earlier in the perturbation approach, 
see Fig^ The nature of this effect is due to the specific 
kx and ky dependence of matrix elements (111(1 . The el- 
ements V^^j^^i which are responsible for the degeneracy 
lifting at kx = ±7r/a and ky = vanish at fcj^ = 0. The 
other set of elements y^,^"^ is non-zero at fc^^ = 0, and 
it opens the gaps inside the BZ. The fcy-dependence of 
the energy bands at fc^; = is shown in FiglJb. The 
degeneracy at ky = is lifted at finite ky by mutual 
influence of linear ky terms in ^ and by the matrix 
elements (|llll . At certain conditions the anticrossing of 
the dispersion curves from different bands 0j may take 




FigOl Energy spectrum at a = 5 • 10^^^ eVm, Vb = 
1.7 meV shown (a) as a function of kx at ky = and 
(b) as a function of ky at = 0. In the latter case the 
anticrossing takes place at the point A for band 2. 

place. An example of the anticrossing effect can be see 
in Fig|23 near the point A for the band 2. Below we 
shall see that the anticrossing leads to the spin flip in 
the {kx,ky) plane for quantum states near the anticross- 
ing point. 

2. SPIN POLARIZATION 

The control on the spin polarization is crucial for 
practical implementation of spintronics. Below we show 
that in the SO superlattice the standing spin wave with 
the period equal to the superlattice period is formed. 
We discuss the space distribution of the spin density for 
the states with different quasimomenta k and the distri- 
bution of spin expectation values in the Brillouin zone. 
The latter describes the mean value of spin polarization 
for the electrons travelling in different directions. We 
have calculated the spin density 

Stic{x,y) = {ip]^)^aiipi^ (12) 

for a quantum state V'k in a given band and after the 
space integration we obtained the vector field of 2D spin 
expectation values (cr2,(k), (jy (k)) in the Brillouin zone: 

a,(k) - i^Pk I <T, I V'k). (13) 

In Figl^lwe show the calculated distribution of {ax,cry) 
for two lowest bands 1 and 2 shown in FigI21 One can 
see that the spin polarization is qualitatively modified 
by periodic potential. The uniform curl distribution of 
spins which typical for the 2DEG with SO interaction 
and without the periodic potential is conserved only 
near the BZ center. It can be seen in Fig. |2K that the 
curl distribution is destroyed at the borders kx — ±7r/a 
of the BZ. The principal difference is that at kx = ±7r/a 
the spins are polarized along x axis and ay — and a 
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^k^ (a) Band 1 








(D ) Band Z 
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Figini Spin polarization shown in one quarter of the 
BZ for of the lowest band 1 (a) and the next band 2 
(b). In the latter case the anticrossing takes place at 
point A, leading to the spin flip in the {kx,ky) plane. 
The parameters are the same as in Fig|H 



new type of singularity appears at = ±7r/a, ky = 0. 
More complicated picture shown in Fig|3|3 takes place 
for the spin polarization in the next energy band 2 
shown in FigI21 We see that at the BZ center the curl 
topology of spin polarization (with reversed angular ve- 
locity since A = 1 for the upper Rashba band) is again 
unchanged. The greatest changes from the uniform curl 
distribution also happen near the borders k^ = ±7r/a 
of the BZ where a new curl has emerged. Another im- 
portant feature of this spin distribution is the spin flip 
at the point A for band 2 shown on the axis kx — in 
the BZ. One can easily establish this point as the anti- 
crossing point for the energy spectrum e = e{kx = 0, ky) 
shown in Figl^b- It should be mentioned that such ef- 
fect was studied previously in quantum wires with SO 
interaction Besides, one can mention in Fig|31that 
at certain points in the BZ the spin average values are 
decreased. The nature of such behavior is the mutual 
influence of SO terms and periodic potential and will be 
discussed below. 

Another effect which can be promising for the ap- 
plication of SO superlattices can be found when an ex- 
ternal electric field in the y direction is applied. Under 
these conditions the spin distribution in the BZ is shifted 
homogeneous along k^ which can produce the magne- 
tization along X direction. It is clear that the sign of 
this magnetization changes with the direction of the Ey 
component of the external electric field and thus can be 
easily controlled. 



3. SPINOR SYMMETRY AND SELECTION 
RULES 

31 Spin parity 

The symmetry of the confinement potential V{x) = 
V{~x) leads to the existence of an additional quantum 
number called the spin parity [1]. Namely, the Hamil- 
tonian commutes with the spin parity operator 

Sx = P.(7cc (14) 

where Px is the inversion operator of the x coordinate, 
Pxf{x) = f{—x). We shall study the cc-dependence of 
the spinor components ip^''^{x) of quantum state l^. 
Below we demonstrate that the Bloch spinors in our 
problem taken at the points k = (fcj, — 0, ±T:/a, ky) can 
be labeled by a certain quantum number s which is the 
spin parity. In this case the spinor components ip^''^{x) 
of © satisfy the following symmetry relations: 

V'k^(a;) = sVk^(-a;), (15) 

where s = ±1 is the spin parity. 

In FigEJjjC we show the x-dependence of the spinor 
components ip^''^{x) in a superlattice elementary cell to- 
gether with the superlattice potential V{x) shown in 
FigS^. One can see that the spinor components ■4!^''^{x) 
are real at kx = ±7r/a. The states are taken at opposite 
ky momentum components and for two neighboring en- 
ergy bands 1 and 2. By comparing Fig0{a) and0Jb), 
01c) one can see that the space symmetry of a partic- 
ular wavefunction on the superlattice period does not 
follow the symmetry of the superlattice potential V{x). 
The reason is the SO coupling which produces the space 
shift of the components ip^''^{x) of the spinor 10. Such 
effect was observed previously for the case of quasi one- 
dimcnsional SO systems ^ . Another important feature 
which can be seen in Figgis that the spin parity changes 
both from the lowest band 1 to the higher band 2 and it 
changes also under the reflection ky —ky. In the next 
Section we shall see how these properties influence on 
the optical selection rules for the interband transitions. 

The symmetry relation H15|l leads to the following 
property of the spin density distribution (|12|1 : 

Sy^zix) = -Sy^z{-X). (16) 

Since the spinor components il^^''^{x) are real at kx = 
±7r/a, the Sy component of spin density vanish iden- 
tically. So, we calculate the space distribution of the 
other components, Sx{x) and Sz{x) taken for the point 
kx = T^/d, ky — iO.OSvr/a in the Brillouin zone of the 
lowest energy band (the corresponding spin polariza- 
tion which is the space integral of Sx,y is shown above 
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(a) 



V(r() = l-'g COS 27rr/a 




(b) Band 1 



■ = 0.05jr/a S=1 
=-0.05jr/fl S= -1 
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(c) Band 2 



-0.057r/a S=1 
= 0.05jr/a s= -1 
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FiglH (a) supelattice potential V(x) and (b), (c) 
space dependence of the spinor components il>^'^{x) in 
a superlattice elementary cell at ky — 0.05±7r/a shown 
(b) for the lowest band 1 and (c) for the next band 2. 
The spin parity s changes from band 1 to band 2 and 



when ky 



ky 



(a) = 

kj, =-0.05?r/fl 



S, ix) 




(b) = 0.05- SCx) = iS/^), S,(x)) 







(C) =-0.05- SCx) = C^:,Cx),;:/x)) 



FigEl (a) spin density distribution Sx{x), Sz{x) and 
(b), (c) the a;-dependence of the two-dimensional spin 
vector {Sx{x), Sz{x)) shown for the lowest band 1 (see 
Figl^J at kx — Tr/a, ky = ±0.057r/a. The spin com- 
ponents switch their signs when the ky momentum pro- 
jection is reflected due to the change of the spin parity 
quantum number. Inside the superlattice elementary 
cell one can observe in (b) and (c) the spin standing 
wave with the spin wriggling in opposite directions for 
the opposite ky components of the electron momentum. 



in Fig|2ta)). It should be stressed that in SO super- 
lattice the Sz component transverse to the superlattice 
plane appears while it is absent in pure Rashba 2DEG. 
The results for Sx{x), Sz{x), and the a:-dependence 
of the two-dimensional vector {Sx{x), Sz{x)) are shown 
in FigO One can see that the average values of the 
spin density across the superlattice cell shown in FigO 
are close to zero which can be seen also in Fig|31in the 
vicinity of the point = ir/a, ky = ±0.057r/a where 
no clear indication of an arrow is seen. In Fig|Sl it is 



6 



V. Ya. Dcmikhovskii, D. V. Khomitsky 



clear also that the spin components change their signs 
when the ky momentum projection is changed, thanks 
to the change of the spin parity quantum number ob- 
served above. In more details this effect can be ex- 
plained with the help of the operator composition for 
the reflection in the reciprocal space. Namely, the trans- 
formation Ky = ky ^ —ky ai the Brillouin zone bound- 
ary kx = ±7r/a is equivalent to the inversion operator 
iiT = k — > — k since the lines kx = ±7r/a are topolog- 
ically identical and one can write KyKx = K where 
Kx = kx ^ —kx and Kx = 1 at the Brillouin zone bor- 
der kx = ±7r/a. Thus, the reflection of ky here produces 
the same effect as the inversion k — > — k, the latter giv- 
ing us the spin flip. As it was mentioned earlier, the 
spin parity and the corresponding symmetry relation 
for it at ky ^ —ky should exist also at kx = 0. Indeed, 
at fca; = the operation ky —ky is also identical to 
k — > — k thus producing the change s ^ — s and making 
a spin flip. 

By looking onto Fig|SIb,c), one can call the spin 
density distribution inside the superlattice elementary 
cell as a standing spin wave. From Fig|3^b,c) one can 
see that the wriggling takes place with the opposite di- 
rection for the opposite ky components of the electron 
momenta in accordance with the symmetry properties 
discussed above. 

32 Selection rules 

The analysis of the optical properties of the SO su- 
perlattices requires the knowledge of corresponding se- 
lection rules. We discuss here the selection rules for the 
operator of the dipole momentum x. The symmetry 
analysis of the wavefunction carried out above will help 
us to define the selection rules for the certain points in 
the BZ kx — 0, zLir/a which are often correspond to the 
energy minima and maxima. By considering the sym- 
metry properties (|15|) of the spinor ^ one can see that 
the matrix element 

^™m'fe = ii^msk I X I lP^,,,k) (17) 

for the direct transitions between the bands m and m' 
having the spin parity s and s' is equal to 

A^^l'fe - (l-ss')AWfe. (18) 

Here Mmm'k is an integral independent of spin parity 
and in general it is non-zero. So, the direct optical 
transitions are allowed only between the states of the 
opposite spin parity. For example, for our problem the 
states in two lowest neighboring spin-split bands 1 and 



2 discussed above satisfy to this condition and the direct 
transitions between them are thus possible. 

In conclusion we would like to mention also that in 
the presence of and external constant electric field par- 
allel to the superlattice direction one can expect a si- 
multaneous appearance of wriggling electron trajecto- 
ries [zitterhewegung ^21) and Bloch oscillations. Such 
situation may lead to non-trivial dynamics and trans- 
port of charged particles, and we plan to investigate 
these problems in a separate paper. 
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